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0\ . 

^\ . We examine global dynamics and bifurcations occurring in a truncated model of a stellar mean 

' field dynamo. This model has symmetry-forced invariant subspaces for the dynamics and we find 

examples of transient type I intermittency and blowout bifurcations to transient on-off intermittency, 
■ involving laminar phases in the invariant submanifold. In particular, our model provides examples 

' of blowout bifurcations that occur on varying a non- normal parameter; that is, the parameter varies 

the dynamics within the invariant subspace at the same time as the dynamics normal to it. As a 
consequence of this we find that the Lyapunov exponents do not vary smoothly and the blowout 
I I bifurcation occurs over a range of parameter values rather than a point in the parameter space. 

> ■ 

i> : 

O . I- INTRODUCTION 

o ■ 

. There has recently been a great deal of work on deterministic dynamical systems possessing invariant submanifolds, 
— . motivated by symmetric systems, and in particular some coupled systems. Such systems have been shov^rn to be 
^\ , capable of producing a range of interesting and robust dynamical modes of behaviour, such as riddled basins and 
on-off intermittency shown to be related by the blowout bifurcation as well as transient on-off intermittency 
We denote the manifold on which such a system is defined by M and the corresponding invariant submanifold 
by Ml. An attractor is said to have a riddled basin if every open set intersecting the basin also intersects the basin 
of another attractor in a set of positive measure; such basins can arise robustly for Milnor attractors contained in 
Mi. 

On-off intermittency to a state in M/ is characterised by dynamics that comes arbitrarily close to the state in Mi but 
^ , that also has intermittent large deviations from M/. Transient on-off intermittency is a transient dynamics exhibiting 
characteristics of on-off intermittent behaviour. Namely, it has an average distribution of laminar phases that satisfies 
a power law with exponent —3/2 [Q. All these phenomena arise as different aspects of blowout bifurcations, where a 
maximum normal Lyapunov exponent of an attractor in M/ passes through zero and thus causes a loss of stability. To 
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?H \ date, most mathematical understanding of such systems is limited to cases where there are a number of simplifying 
. . 1 assumptions: 

(HI) The control parameters are assumed to be normal in the sense that the dynamics of the invariant sub- 
manifold is independent of the bifurcation parameter. Such parameters preserve the dynamics on the invariant 
submanifold, but allow it to vary in the rest of the phase space; 

(H2) The attractor which becomes transversely unstable in the Mi and hence causes the blowout bifurcation is usually 
assumed to be chaotic. 

(Note however that j?) find blowout type bifurcations from quasiperiodic dynamics on T^.) 

As many physical systems of interest are unlikely to possess normal parameters (a notable exception to this being 
some coupled systems), we expect that (HI) is not usually applicable. This is particularly expected to be true in 
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truncations of partial differential equations that arise in fluid and dynamo equations (as well as astrophysical models 
in general). Thus, by studying the behaviour of a specific example where (HI) does not hold, we hope we throw some 
light on the dynamics of general systems of this type. 

To this end we examine a system where (HI) does not hold, i.e. where the control parameter varies the dynamics 
within the invariant subspace as well as that normal to it. We see that this has the effect of 'spreading' the blowout 
bifurcation out over an interval of parameter values due to the existence of periodic windows where (H2) does not 
hold; however, we conjecture that there is a positive measure subset of parameters on which the blowout resembles 
the case for normal parameters. 

The model we describe in Section H arises as the truncation of a stellar axisymmetric mean field dynamo model 
where the natural control parameters are not normal. There is also more than one invariant manifold forced by the 
spatial symmetries of the system, although this does not appear to affect the behaviour in the cases examined, in the 
sense that it is only one of these manifolds, namely the antisymmetric one, which seems to dominate the attracting 
dynamics. 

In Section III we discuss numerical results from simulations of the model; we discuss the basic bifurcational behaviour 
in the full system and the antisymmetric subspace before discussing examples of type I intermittency, transient on-off 
intermittency and non-normal parameter blowout bifurcation in the model. In the final Section IV the implications 
of the results are discussed for more general systems of this type as well as for the dynamo problem. 



II. MODEL 



A. Galerkin model for the mean field dynamo 

The dynamo model considered here is the so called ail mean field dynamo model, with a dynamic a-effect, given 
by Schmalz & Stix (see also Covas et al. for details). We assume a spherical axisymmetrical configuration with 
one spatial dimension x (corresponding to the latitude and measured in terms of the stellar radius R) for which the 
magnetic field takes the form 

1 dA, 



where A is the (/)-component (latitudinal) of the magnetic vector potential and B is the (/)-component of B. The 
model is given by the mean field induction equation 

dB 

— ^yx{vxB + aB- r/fV x B), (1) 

where B is the mean magnetic field, v is the mean velocity, rjt is the turbulent magnetic diffusitivity and a is the 
coefRcient of the a-effect In addition the a-effect, which is important in maintaining the dynamo action by 

relating the mean electrical current arising in helical turbulence to the mean magnetic field, is assumed to be dynamic 
and expressible in the form a = aocosx — aM{t), where ao is a constant and um is its dynamic part satisfying the 
equation 

^ = "*^ + ^^-^' ^2) 

where Q is a physical constant, J is the electrical current and Vt is the turbulent diffusivity. These assumptions allow 
Eq. (n|) to be split into the following two equations: 



dA^ ^ d^A^ 
dt i?2 dx^ 
dB^ rjt d^B^ luq 9 A, 



dt i?2 Q^2 • R dx ■ 
Expressing these equations in a non-dimensional form, relabelling the new variables to 

(^0, B^, um) => [A, B, C), 

and using a spectral expansion of the form 



aB^, (3) 

(4) 
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TV 

A — An(t) sin nx, 

n=l 
N 

B = Y,Bn{t) Sin 

n=l 
N 

C7 = ^C„(t)sin 



nx, 

n=l 

where N determines the truncation order, reduces equations (||)-(||) to a set of ordinary differential equations, the 
dimension of which depends on the truncation order N . We consider the fuU system given in terms of the variables 
Bn, Cm n = 1, • • • , iV in the form 

8 A D 

^ = -n'A,, + j{Bn-i + (5) 

N N 

+ ^ ^F(n,m,0-B™Cz, 

m=l 1=1 

N 

-n^Bn+^G{n,m)A^, (6) 



dBn 

dt 
dt 

where 



ni— 1 

N N 



-vn^Cn -^Y^ H{n, m, l)A,nBi, (7) 



m=l 1=1 



F(n, m, I) — 

8nml 



7r (n+m+^ ) (n+m—/){n—m+/) (n — rn—T) ' 

if n + m + / is odd and F(n, m, I) — otherwise, 

H{n, m, I) — 

4 nml{-n^+3m^+l^) 



if n + m + Z is odd and H{n, to, Z) = otherwise and 

4:71772 

G{n,m)^—-^ JT' (8) 

7r(7i^ — m^j 

if 71 + TO is odd and G'(77, to) = otherwise. In these equations the control parameters are the so called dynamo 
number D (which is proportional to the square of the angular velocity gradient and to the square of the turnover time 
of the turbulent convection eddies) and the diffusivity ratio v = 



B. Invariant subspaces for the model with A'^ = 4 

Covas et al. [|j confined themselves to the study of models that are antisymmetric with respect to the equator and 
found that the minimum truncation order N for which a similar asymptotic behaviour existed was iV = 4. In this 
case, the equations have a twelve dimensional phase space and are symmetric under the four-element Abelian group 
that comprises the identity /, the reversal transformation 

Ar,{t) ^ -An{t), B„{t) ^ -Bn{t), Cn{t) Cn{t), 

the antisymmetric (or dipolar) transformations 

A„(i) ^ (-1)"+M„(t), Br,{t) ^ (-l)"i?„(t), C„(t) ^ (-l)"C„(i) 
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and the symmetric (or quadrupolar) transformations 

The trivial solution, given by An — Bn = Cn ~ 0, is the only one which possesses both the dipolar and quadrupolar 
symmetries while symmetric solutions come in pairs and asymmetric solutions come in quadruples. 

The antisymmetric part of these equations, which is obtained by putting Bi — Ci — A2 — B^ — — A^^ = 0, was 
studied in . We refer to this dynamically invariant subspace 

= 0,0, 0,^2, 6*2,^3,0,0, 0,^4,^4)} 

as the antisymmetric subspace. There is also a six dimensional symmetric invariant subspace 

Ms = {(0, Si, 0, A2, 0, C2, 0, B3, 0, A4, 0, C4)} 

although as we will see, the attractors are typically not contained within Ms- Throughout the paper we refer to the 
full system as the 12 - dimensional system. 

III. DYNAMICAL BEHAVIOUR 

The system considered here has a two dimensional parameter space {D, v), neither of which is normal for the system 
restricted to Ma, as can be seen from equations (||)-(0). We confine v to the range [0, 1] on physical grounds, as 
otherwise there will be no dynamo action. Previous studies of these models have taken 1/ — 0.5. Here we shall consider 
two distinct cases of v given by 0.5 and 0.47 in this range and in each case allow D to vary. To study the dynamics 
of this system, we start by looking at the dynamics on the antisymmetric invariant submanifold Ma and then look 
at how this changes as the full system is switched on. 

A. Basic bifurcation behaviour 

To begin with, we consider the case of = 0.5 and as a first step make a coarse study of the dynamics confined to 
Ma as well as the full (12-dimensional) system by considering the averaged energy (E cx \B\'^dx) as a function of 
the parameter D. The results of these calculations are summarised in Figs. (|]) and (||) respectively. The figures were 
produced using a fourth order variable step size Runge-Kutta method to integrate a number of randomly selected 
initial conditions forward in time, and so get a selection of the possible attractors. After a time when transients were 
deemed to have decayed (which we took to be 1000 time units) we averaged the energy over a much longer time series 
i.e. 10000 time units. We have verified the following results using the continuation package AUT097 Q. 

For small D {D < 98.67) all attracting dynamics of the full 12-dimensional system is confined to the 6-dimcnsional 
antisymmetric invariant submanifold Ma- The details of bifurcations are depicted in Figs. (^ and (^). As can be seen 
from these figures, as D is increased, the fixed point at the origin (the trivial solution for both systems) bifurcates at 
D = 12.57 to two fixed points, which are symmetric with respect to An — > —An, Bn — > —Bn, Cn — > Cn- A subcritical 
pitchfork-type bifurcation stabilises the origin from D — 77.25 up to a supercritical Hopf bifurcation at D — 98.67. 

For larger values of D the two systems will evolve differently. In AIa the supercritical Hopf bifurcation creates 
a stable periodic solution while the 12-dimensional system evolves on a branch of stable periodic solution with 
quadrupolar symmetry that becomes unstable through a torus bifurcation at D = 118.76. In addition to this, at 
D — 118.23, the origin undergoes Hopf bifurcation creating a branch of dipolar periodic solutions that are stabilised 
by a torus bifurcation at D = 127.78, as shown in the inset of Fig. (|^). Also shown in the inset is the appearance 
through a saddle-node bifurcation of a branch of periodic solutions in Ms at D = 132.72 that are stable up to 
D = 135.42. Here they lose stability at a torus bifurcation. All of these branches of solutions pass through the region 
from D — 150 to D — 175 without bifurcation; however they do not provide a complete picture of all that happens in 
this region. 

There is a pair of periodic orbits in Ma that are created at a saddle-node bifurcation at D — 170.25. These are 
stabilised by a pitchfork bifurcation of periodic orbits at D = 171.10 creating a branch of stable periodic orbits with 
no symmetry continuing down to D = 171.003; these are important for the intermittent dynamics discussed in the 
next subsection (Fig. ||). For larger D, the periodic orbit undergoes a sequence of bifurcations preserving the dipolar 
symmetry (i.e. creating branches that remain within Ma) up until D ~ 111 -7b where a branch of stable periodic orbits 
that bifurcate from Ma is created. We conjecture that the saddle node bifurcation creating this series of periodic 
orbits is associated with breakdown of a quasiperiodic flow on a two-torus created at one of the torus bifurcations but 
have not been able to check this. 
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B. Intermittent dynamics 



We investigated two cases; v — 0.5 and v = 0.47. Note that for physically meaningful results we require v < 0.5. 
We have examined the transverse stability of attractors in Ma by calculating the corresponding transverse Lyapunov 
exponent (At). Fig. ^ summarises the results of calculations of the At for the periodic orbit which is created at 
D = 170.25, as shown in the inset in Fig. (|l|). The important feature in this behaviour of At is the presence of the 
two crossings through zero. The transverse stability of the other attractors do not change in this range of D. 

1. The case of v — 0.5 

For the case v — 0.5 we examined the behaviour of the full system over a parameter interval in the neighbourhood 
of -D = 171. The loss of stability of the periodic orbit described in the previous section, docs not induce on-off 
intermittency as we first suspected. Instead, the bifurcation at D — 171.10 is a pitchfork bifurcation that creates an 
asymmetric periodic orbit which survives up to Z? = 171.003 and is then destroyed by a collision with an unstable 
orbit in a saddle node bifurcation, which is shown in Fig. ^. 

For D < 171.003, we have a transient type-I intermittency, as can be seen in Fig. ||. We also calculated the scalings 
of the transient times and average times between the bursts (as shown in Fig. g and Fig. 0) and found them to be 
in good agreement with the known —1/2 scaling. The behaviour between bursts shown in Fig. ^, resembles the 12D 
periodic orbit discussed above, except that the amplitude of the symmetric part of the oscillations between the bursts 
grow slowly towards the bursts and return, after the bursts, close to the invariant submanifold. The intermittent 
behaviour is transient, in the sense that the orbit returns to the fixed point in the invariant submanifold. The spectrum 
of Lyapunov exponents for these transient forms of intermittency is in the form (+,0, — ...), until the transient dies 
out, becoming attracted to a stable fixed point (note: there are two such fixed points, located symmetrically with 
respect to the A and B variables). 

On the negative side of the crossing of the transverse Lyapunov exponent, shown in Fig. ^, we observe a basin 
boundary for the full 12-dimensional system with a dimension close to that of the phase space. This is shown in Fig. 
P which demonstrates which asymptotic attractor on the invariant submanifold the initial conditions get attracted 
to. Both basins seem to be made up of open sets (supported by the fact that calculations indicate an integer box 
counting dimension). This conclusion is further supported by the calculation of the exterior dimension |l3| shown in 
Fig. |. 

It is also interesting to note, from both physical and mathematical points of view, that even if the full (12- 
dimensional) system does have any new attractors, nevertheless it will in general alter the relative size of the basins 
of attraction; most initial points seem to get attracted to only the fixed points or one of the cycles, as opposed to 
initial conditions starting in the invariant sub-manifold. 

2. The case of u = 0.47 

By making v slightly smaller than 0.5 we were able to change the order of the bifurcation to chaos in the invariant 
submanifold relative to the loss of transverse stability. To study the behaviour of the system with v — 0.47, we looked 
at the parameter region in which the system with = 0.5 had a chaotic attractor (as depicted in Fig. For our 
calculations, we chose our initial conditions to lie in the basin of the chaotic attractors for = 0.5 system. We then 
studied the evolution of the system for v = 0.47 by changing the control parameter and taking the initial conditions 
at each step to lie in the basin of the attractor for the previous parameter value. The results of these calculations for 
the At and the largest Lyapunov exponents of the full 12-dimensional system are given in Fig. |l^ and Fig. |l|. As 
can be seen, the chaotic behaviour is now interspersed with periodic windows. Within these windows the periodic 
solutions co-exist with chaotic repellers. 

Another crossing of the transverse Lyapunov exponent, from negative to positive, shown in Fig. ^, occurs at 
D — 178.71, and for 177.10 < D < 178.71 there are stable periodic orbits restricted to the invariant submanifold. 
This crossing is also related to a bifurcation of periodic orbits and therefore there is no indication of on-off intermittency 
or riddled basins. The periodic orbit on the invariant submanifold becomes chaotic just after D = 178.76, not close 
enough to the transverse stability bifurcation at I? = 177.71 to induce on-off intermittency. This suggests that there 
is likely to be a blowout at nearby parameters in the two parameter space, as we discuss in the next section. To 
substantiate this, we calculated the scaling of the probability distribution of the off phases (corresponding to when the 
distance to the invariant submanifold is less than 10"'^) as a function of their length. As can be seen in Fig. O, the 
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scaling agrees with the power law behaviour proposed in [g| with an index of —3/2. At this parameter value transient 
on-ofF trajectories appear to be induced by a chaotic invariant set that is a repeller within the invariant submanifold. 



C. Generic behaviour for non-normal parameters 

There are a number of interesting dynamical phenomena that occur here, that are related to the fact that the 
system parameters are not normal. 

(I) The chaotic behaviour in the invariant submanifold appears to be of the non-uniformly hyperbolic variety, and 
in particular the chaotic attractors are not structurally stable; the y a re destroyed by arbitrarily small perturbations. 
Notwithstanding this, we find numerical evidence (Fig. ^and Fig. |ll|) that there is a family of chaotic attractors with 
similar properties defined on a subset of parameter space with positive (Lebesgue) measure but open (even dense) 
complement. This is what is found, for example, in the logistic map In the open dense complement we expect to 
see periodic windows and many bifurcations, for example period doubling cascades, which we have found numerically. 
This is consistent with the conjecture of Barreto et al. (l6j on noting that the attractors here have only one positive 
Lyapunov exponent. In this parameter region the system may be said to be fragile in the sense that arbitrarily 
small changes in the control parameter D can force a chaotic attractor to be replaced by a nearly attracting periodic 
orbit. 

(II) In the light of (I), there is no reason why there should be a unique parameter value Dc at which blowout occurs. 
In particular, the attractor in Ma varies discontinuously, and its tangential and normal Lyapunov exponents vary 
discontinuously with D except within the periodic windows. This explains the presence of smooth segments in the 
curves of Fig. [llj within regions where the attracting dynamics is periodic. 

(III) In this system the passage of At through zero is fairly simple. Firstly, the value of At can be bracketed 
between upper and lower bounds that also pass through zero. This is presumably due to the fact that certain periodic 
orbits in the attractor will typically maximise and minimise transverse Lyapunov exponents [^5|. Secondly, there is 
evidence that there is a positive measure Cantor set S in parameter space where chaotic behaviour exists. On S the 
Lyapunov exponents are continuous in the sense that there is a continuous function A' of parameter that is equal to 
At on S and passes through zero at about D = 177.75 (see Fig. |lO| and also Fig. Even for parameter values not 
in iS we can get transient on-off intermittent behaviour (see Fig. [L2D . 

On the basis of our results we conjecture that properties (I)- (II) are typical behaviour at blowout on varying a 
non-normal parameter and (III) is a typical simple scenario of how this can occur. 

For the case v = 0.47 we note that the transverse Lyapunov exponent (of the chaotic invariant set that attracts 
within Ma) becomes positive and causes the appearance of transient on-off intermittency. 



IV. DISCUSSION 



We have studied global dynamics and bifurcations occurring in a 12-dimensional truncation of a stellar mean field 
dynamo model which possesses two six dimensional invariant submanifolds corresponding to dipolar and quadrupolar 
symmetries respectively. An essential feature of this model is that its control parameters are non-normal, allowing 
the dynamics to vary both within the invariant submanifolds as well as in the directions normal to them. 

Depending upon the region of the parameter space considered, we find a diverse set of dynamical modes of behaviour, 
including different forms of intermittency. In addition to transient type I intermittency, we find transient on-ofF 
intermittency induced by blowout bifurcations. In the parameter range where we observe the latter behaviour, the 
invariant submanifold possesses a family of chaotic attractors on a subset of parameter space with positive (Lebesgue) 
measure but open (even dense) complement. On the basis of our numerical calculations these attractors seem to be 
structurally unstable, which is consistent with the conjecture of Barreto et al. |l^]. We also find that as a consequence 
of the non-normality of the control parameters the blowout bifurcation seems to occur over an interval rather than a 
point in the parameter space. 

These results can be of potential significance for the dynamical behaviour of systems with non-normal parameters. 
Given the fact that the model considered here was derived directly from dynamo equations, the forms of intermittency 
found here can also be of potential importance in understanding the mechanism of production of the so-called grand or 
Maunder-type minima in solar and stellar activity, during which the amplitudes of stellar cycles are greatly diminished 
|p8| , p^ . We do not however wish to imply that the forms of intermittency responsible for such stellar behaviour are 
necessarily transient. 
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FIG. 1. Diagram showing attractors for a random selection of initial conditions within Ma for v = 0.5. TFP stands for 
trivial fixed point, FP for non-trivial fixed point, PA for antisymmetric periodic orbit and CA for antisymmetric chaotic orbit. 
Continuation using DSTOOL shows that the break in the FP attractor uppermost in this diagram is just a feature of the choice 
of initial conditions; in fact it continues to be attracting over the whole range of D. The inset shows coexisting chaotic and 
periodic attractors over a range of D. 

FIG. 2. Diagram showing attractors for randomly chosen initial conditions in the full phase space for v = 0.5. PS stands for 
symmetric periodic orbit (i.e. in Afg), PM for a periodic orbit neither in Ma nor Ms and QPA for antisymmetric quasiperiodic 
orbit. The notation of TFP, FP and PA are as in Fig. |l| Observe the existence of intermittent behaviour over a range of D. 

FIG. 3. The leading (i.e. most positive) transverse Lyapunov exponent and the two leading Lyapunov exponents for the 
attractor of a particular initial condition for the system on Ma at v = 0.5. The attractor is a periodic orbit from D = 170.25 
up to -D = 178.76, although it undergoes a number of period doublings in this range to give a chaotic attractor for D > 178.76. 
The computed orbit is transversely stable in the range 171.10 < D < 178.71. By reducing u to 0.47 we can change order of the 
loss of transverse stability and the breakdown to chaos in Ma- 

FIG. 4. Continuation of a periodic orbit showing breakdown to type I intermittency at = 0.5. The abscissa shows a 
symmetric component of a branch of PM periodic orbits created at a pitchfork of a PA periodic orbit. This is destroyed at a 
saddle-node bifurcation giving rise to type I intermittency at D < 170.003. 

FIG. 5. Time series showing a component transverse to Ma for transient type I intermittency series for u = 0.5 at D = 170. 
Observe the long but irregular periods of lingering near a small amplitude periodic orbit interspersed by large fluctuations. 
After a long time, the trajectory is asymptotic to a stable fixed point. 



FIG. 6. Scaling of transient time of the transient type I intermittency for v = 0.5 against the unfolding parameter D — D^. 



FIG. 7. Scaling of transient burst time of the transient type I intermittency for v = 0.5. 

FIG. 8. Two dimensional slice through phase-space obtained by setting all components zero except for Ai and Bi. The 
basins of attraction of the fixed point (black) and the periodic orbit (white) that are coexisting attractors at D = 171.12 and 
V = 0.5. 

FIG. 9. Approximation of the exterior dimension of the basin shown in black in Fig. ^ This is very close to the 
dimension of the slice through phase-space indicating that the basin boundary is highly convoluted, even though it is not 
riddled. 

FIG. 10. Largest transverse Lyapunov exponent (At) and the two leading Lyapunov exponents (Ai and A2) for a family of 
attractors on the antisymmetric invariant submanifold for u — 0.47. Note the existence of periodic windows and general trend 
of At through zero indicate a blowout bifurcation near D « 177.75. The lack of smoothness of these curves is indicative of the 
fact that D is not a normal parameter. 

FIG. 11. Amplification of the transverse Lyapunov exponent and the two leading Lyapunov exponents for the antisymmetric 
subset of equations for f — 0.47. This shows the existence of a "window" in parameter space where the attractor within the 
invariant subspace is periodic and transversely repelling. 

FIG. 12. Time series showing transient on-off intermittency for v = 0.47 and D = 177.70. The on-off intermittent behaviour 
is induced by a chaotic repeller which is present within the periodic windows. After a long transient, the trajectory is asymptotic 
to a stable fixed point within Ma (not shown). 
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FIG. 13. Scaling of the laminar phases over an on-off transient orbit segment for i/ = 0.47 and D — 177.70. The —3/2 
scaling is evidence of an on-off intermittent state. 



FIG. 14. Average of the variable B2 measuring average distance from Ma over an on-off transient orbit segment that 
eventually ends at a fixed point (for f — 0.47). The discontinuous nature of this presumably reflects the discontinuous change 
in At illustrated in Figure |l^. 
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